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Abstract 

The number of neutrino induced upward going muons from a single gamma 
ray burst (GRB) expected to be detected by the proposed kilometer scale Ice- 
Cube detector at the South Pole location has been calculated. The effects of 
the Lorentz factor, total energy of the GRB emitted in neutrinos and its dis- 
tance from the observer (redshift) on the number of neutrino events from that 
GRB have been examined. The present investigation reveals that there is pos- 
sibility of exploring physical processes of the early Universe with the proposed 
kilometer scale IceCube neutrino telescope. 
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1 Introduction 



We consider the relativistically expanding fireball model of gamma ray bursts (GRB) 
|T[] and calculate the number of secondary muons from high energy neutrinos of a 
GRB expected to be detected in a kilometer scale neutrino telescope. In a relativistic 
expanding fireball protons may be accelerated by Fermi mechanism to energies as high 
as lO'^^eV. The accelerated protons interact with photons of energies of the order of 
MeV to produce pions. From protons both 7r"'"'s and vr^'s are produced sharing the 
proton energies roughly in equal parts. The tt^'s would lead to generation of high 
energy neutrinos and from vr'^'s high energy photons would be produced. 

The decay process for neutrino production from 7r"^'s is as follows tt"*" ^ yU"*" + 
z/^ + Vf. + + The Lorentz factor F of a GRB has a crucial role in the 

neutrino production mechanism of the GRB. The efficiency for producing pions in 
p — 7 collisions in the fireball varies as F~^ and the break energy for neutrino energy 
spectra varies as F^. Although F grows linearly with the radius in a relativistic fireball 
model, F saturates at a value of the order of 100 [|l|. 

The fraction of energy lost by protons to pions /,r can be expressed as a function of 
the pair production optical depth r^^. If the optical depth is large there would be 
more pion productions. Also if the optical depth is large then the photons produced 
from 7r° decays can not come out as a burst. In that case the GRB would be opaque 
for photons. The GRBs which are not bright in photons are bright in neutrinos. 

The effects of varying the Lorentz factor F, the total energy emitted or released in 
neutrinos Egrb and redshift of the GRB z on the expected number of neutrino in- 
duced muons from a GRB have been studied in this work. Each gamma ray burst is 
an individual phenomenon in the Universe distinguished by its redshift, total energy 
emitted in neutrinos or photons, duration of the burst and Lorentz factor at a partic- 
ular instant of time . To obtain informations on the physical process responsible for 
a gamma ray burst it is always meaningful to carry out studies on the detectability 
of individual GRBs. 

In Ref. [0] the correlations of the parameters minimum Lorentz factor F, wind 
variability time t^, observed photon spectral break energy E^^j^^y for wind luminosity 

= 10^^ erg are given in FIG. 1. for equipartition parameters sb = 0.01 and = 0.1. 
We have taken points from the contour plots given in FIG.l. of Ref. and used 
those values in our calculations. In the present work we use the wind luminosity 
but by Egrb we mean the total energy released in neutrino emission. The observed 
photon spectral break energy E!^ j^^y is the energy where luminosity per logarithmic 
photon energy interval peaks. 

The fraction of fireball energy that goes into neutrino production is weakly de- 
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pendent on the wind model parameters. There are two reasons for that. The first 
reason is for low values of Lorentz factor F and wind variability time only a small 
fraction of the pion's energy is converted to neutrinos at high proton energy due to 
pion and muon synchrotron losses. The second reason is the observational constraints 
imposed by 7-ray observations imply that the wind model parameters F, wind lumi- 
nosity and wind variability time are correlated. We keep in mind the luminosity 
mentioned in is the wind luminosity and in this work also we use the same nota- 
tion. To avoid confusion we mention that the fraction of proton energy that goes into 
neutrino production does not appear explicitly in the expression of total number of 
neutrinos emitted from a GRB in our procedure of calculation. The reason behind 
this is we have used Egrb which is the total energy released in neutrinos to calculate 
the number of neutrinos emitted from a GRB. 

From one proton two muon neutrinos are expected to be generated if both the pion 
and the muon decay. The neutrinos carry about 5% of the proton energy, it is given 
in 0. So the total energy Egrb released in neutrinos of a GRB is about 10% of 
the original total fireball proton energy. We have used the total energy emitted in 
neutrinos Egrb to calculate the neutrino spectrum from a GRB. So the energy loss by 
protons in pion productions as well as the pion and muon energy losses which are the 
intermediate processes for neutrino productions from protons have been already taken 
into account in our procedure of calculations. The small value of Egrb compared to 
the original total fireball energy accounts for the intermediate energy loss processes 
in production of neutrinos from protons in a GRB. If we use Egrb = lO^^erg in our 
calculation that means the total fireball proton energy was lO^^'erg assuming 10% of 
the fireball energy has been released in neutrino emission from the GRB. The wind 
durations are of the order of 10 seconds. 

The neutrino spectrum from a GRB on Earth depends on the total energy emitted 
in neutrinos Eqrb, distance of the GRB from the observer z and the neutrino break 
energy E^ which is again dependent on the Lorentz factor F, observed photon spectral 
break energy E^ j^j^y, the wind variability time as well as wind luminosity L.^ and 
equipartition parameter e^. 

Earlier average neutrino event rates in a neutrino telescope of km? area per year 
from gamma ray bursts have been obtained considering burst-to-burst fluctuations in 
distance and energy ^ . They have used energy and redshift distribution functions 
of the GRBs modeled according to observations on GRBs. The neutrino event rates 
per year from GRBs have been obtained including the fluctuations in the Lorentz 
factor F. The purpose of the present work is not to calculate the neutrino event rate 
expected per year in a muon detector of km? area. We investigate for what range of 
values of their physical parameters {Egrb, z, F) individual GRBs will be detectable 
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by a kilometer scale neutrino telescope of a given threshold energy. 

In GRBs z/^ + z/^'s are almost 10^ times more abundant than + v^. The primary 
neutrino flux from GRBs contain very small number of z/^ + z/^. Due to neutrino 
oscillations some of the z/^ 's convert to z/^-'s during their propagation from the burst 
location to the detector on Earth . According to SuperKamiokande measurements 
the flux ratio of muon neutrinos and tau neutrinos is F^^j^^ = 0.5. The Ur + i^r 
events would produce different experimental signatures from the z^^ + z/^ events [Q. 
Cosmic tau neutrinos while coming close to the surface of the detector may undergo 
a charged current deep inelastic scattering with nuclei inside or near the detector and 
produce a tau lepton in addition to a hadronic shower. This tau lepton traverses 
a distance, on average proportional to its energy, before it decays back into a tau 
neutrino and a second shower most often induced by decaying hadrons. The second 
shower is expected to carry about twice as much energy as the flrst and such double 
shower signals are commonly referred to as double bangs. The tau leptons produced 
in this way are not expected to have any other relevant interactions as they are losing 
energy very fast and subsequently they decay to muons. Other than double bang 
events there are Ur + i^t events in which muons would be detected from r decays 
{ur ^ r ^ fi). The double bang events are much less compared to the events in 
which tau leptons decay to muons (z/^ ^ r ^ /i). The z/^ + z/^ events are much more 
in number compared to the Ur + Vr events. The relative magnitudes of rates of double 
bang events from z/^'s, events in which r's decay to /x's and the + events per 
year in a kim? area neutrino telescope are given in ||^. The majority of muon events 
in the kim? area muon detector would be from interactions of + z/^ in rock (in case 
of upward going muons) or ice (in case of downward going muons) at the South Pole 
location. 

In this work we consider only the upward going muon events from z/^ + z/^ of 
a gamma ray burst because incase of upward going events the background noise is 
comparatively less. 

Recent observations on GRB afterglow have conflrmed the relativistic flreball model, 
to be more speciflc the internal-external shocks model . If GRB neutrinos could be 
observed by neutrino telescopes it would be a great achievement for understanding the 
distribution in Lorentz factor F, the ultrahigh energy neutrino background, testing 
special relativity and the equivalence principle. 

In section 2 we discuss about the parametrization of + z/^ spectrum from a GRB 
and the procedure used here to obtain the number of neutrino induced muon signals. 
In section 3 we mention about the atmospheric neutrino background which often 
makes the detection of neutrino signals from astrophysical sources difficult in muon 
detectors. In section 4 the results of this work have been discussed. 
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2 Parametrizations of the Neutrino Spectrum of 
a GRB and the Secondary Muons Produced in 
Rock 

2.1 Neutrino Spectrum of a GRB 

A total amount of energy Egrb is released in neutrino emission from a GRB at a 
distance of redshift z from the observer on Earth. 

The Lorentz factor of the GRB is F. The neutrino spectrum is expected to be 
correlated mainly with the observed gamma ray spectrum of the GRB . In analogy 
with observed 7-ray spectrum of a GRB we write down the neutrino spectrum of the 
GRB. 

^ = Axmm(l,E./E^)-l (1) 

Here A is the normalisation constant and E^, is the break energy of the neutrino 
spectrum which is a function of Lorentz factor F, photon spectral break energy E^ n,j^y. 
Since F, E^^j^^y and wind variability time are correlated so E^, is also dependent 
on tv for a given value of wind luminosity and equipartition parameter e^. 



K ^ IQ'^^GeV (2) 

In the above expression F2.5 = F/10^'^. E^rnin and E^^nax are the lower and upper 
cutoff energy of the neutrino spectrum. We assume E^rnin « E^. 

The normalisation constant A is determined from energy considerations 

^ = . wf" ^ (3) 

To make things clear we again mention that all the intermediate energy loss terms for 
neutrino productions from protons are accounted for within Egrb which is a small 
fraction of the initial fireball proton energy. The observed neutrino energy on Earth 
Eyobs and its energy at the source, E^, are related as E^obs = E^/ (1 + 2;). The observed 
cutoff energy of the source is E„rnaxobs = Eymax/ (1+2)- We assume Eyrnax = lO^^Gel^. 
Assuming isotropic emission from the source, the total number of neutrinos from a 
single GRB at redshift z per unit energy at energy E^obs that strike per unit area of 
the Earth is 

dN^o dNy 1 



dE^obs dEy 4:7ir'^{z) 



;i + (4) 
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r{z) is the comoving radial coordinate distance of the source. For a spatially flat 
Universe (which we shall assume to be the case) with Qj>^ + = 1) where and 
Qm are, respectively, the contribution of the cosmological constant and matter to the 
energy density of the Universe in units of the critical energy density, 3Hq/8ttG, the 
radial coordinate distance r{z) is given by 

rw^rm , (5) 



Here c is the speed of light, Hq is the Hubble constant in the present epoch. In 
our calculations, we shall use = 0.7, Vim = 0.3, and Hq = 65 km sec~^ Mpc~^. 
A similar procedure for TeV photon flux calculation on Earth from GRBs has been 
used in [0. 

The orginal muon neutrinos from GRBs are undergoing oscillations during propaga- 
tion and tau neutrinos are produced in this way. The probablity of oscillation from 
to is given by 

Prohyv^ — > z/^) = sin 29sin ^ ) (6) 

Considering the values of neutrino mass difference and mixing as found from 
SuperKamiokande data {sin'^29 ~ l,5m^ ~ lO^^el^^) and the distance d ~ 100 
Mpc between AGN and our galaxy the above expression for probablity averages to 
about half for all relevant neutrino energies to be considered for detection [^]. The 
high-energy muon neutrino flux is suppressed by a costant factor of half independent 
of its energy due to neutrino oscillations. The expression for muon neutrino flux from 
GRBs we have obtained above has been divided by a constant factor of two to include 
the effect of neutrino oscillations. 



2.2 Number of Secondary Muons Produced from Muon Neu- 
trinos of a GRB 

The muon events from GRB neutrinos have been calculated including muon energy 
loss in rock and absorption of neutrinos during their passage through the Earth. 
The number of neutrino induced muons from a GRB of energy Egrbi redshift z and 
Lorentz factor F has been calculated above a muon threshold energy E^h using the 
standard procedure from |^. The total number of secondary muons from the muon 
neutrinos of a GRB below the detector above a threshold energy Eth is given by 

S{>Eth)= dE,obs-TT^Pu{Euobs,> Eth)exp[-atot{E,obs)NAX{9)] (7) 
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In the above expression atot is the total cross section for neutrino absorptions. It 
is the sum of charged- current and neutral-current cross sections DTO). X{6) is the 
matter traversed in gm/cm^ by the neutrinos emitted from the GRB during their 
passage through the Earth at a zenith angle 9. Na = 6.022 x 10^^ cm~^. The proba- 
blity that a muon produced in a charged- current neutrino nucleon interaction arrives 
the underground detector with an energy above the muon threshold energy Eth is 
P^{E^obs,> Eth)- It depends on the neutrino nucleon charged-current cross section 
of muon production and the energy loss rate of muons in rock. We have used the 
neutrino- nucleon interaction cross sections from [0. The rate of muon energy loss 



in rock has been used from |TT| . The expression for probablity of muon production 
from neutrinos is 

Pu{E,obs, > Eth) = / dE ^ ^ / dy' / dx-— (8) 

J Eth (e + E^) Jo Jo dxdy' 

The charged-current differential cross sections can be expressed in power laws in 
neutrino energies. 

BEl,, (9) 



dEyQl)g 

where B is the amplitude and 5 is the differential index of power law. The expression 
of probablity of muon production from neutrinos can be simplified to 

> E.,) = ^ j^^EitL\l - (10) 

d + 1 J Eth (e + E^) Euobs 

The expression for rate of energy loss by muons is 

f^- — ''^^ (11) 

where for rock 

a = [2.033 + 0.077 In E^{GeV)]10-^GeVgm-^cm^ (12) 

and 

p= [2.229 + 0.2lnE^t{GeV)]10-^gm-^cm^ (13) 
from |]lT|. In the expression of probablity e = a/(3 and ^ = l/f3. 



The expression for probablity has been substituted in equation (8) and the number 
of neutrino induced muons above muon threshold energy lOOGeV has been calculated 
for a detector of km 
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Atmospheric Neutrino Induced Muons in a Neu- 
trino Telescope of krn^ Area 



The atmospheric neutrino fluxes from |12| and the neutrino nucleon interaction cross 



sections from |]Tol have been used to calculate the neutrino induced muon fluxes from 
atmospheric neutrinos in a muon detector of km'^ area above a muon threshold energy 
Eth = lOOGeV^. The secondary muon fluxes produced in rock have been calculated 
for cos 6' = —1 and cos^^ = —0.05 where 6 is the zenith angle . The procedure 
for calculation of neutrino induced muons as discussed in section 2 is used in this 
case also including the absorption of neutrinos during their propagation through the 
Earth. The IceCube detector would be placed at a depth of 2 km below the surface 
of Earth. We have used the effective matter traversed by neutrinos inside the Earth 
when calculating the neutrino events for cos6 = —0.05. 



4 Results and Discussions 

In FIG. 1. of Ref. 0] contour plots of observed photon spectral break energy are given 
for Lu, = lO^^erg/sec and equipartition parameter = 0.1, = 0.01. From that 
figure for a value of minimum Lorentz factor F we can find the allowed value of E!^ j^j^y 
and the corresponding value of wind variability time ty sec for Ly^ = lO^^erg/sec and 
equipartition parameter = 0.01. For one value of F there can be many allowed 
values E!^ MeV '^^^ values of minimum Lorentz factor F, wind variability time and 
observed photon spectral break energy E^^j^^y used in our work for wind luminosity 
Lyj = lO^^erg/sec and equipartition parameter = 0.01 have been listed in TABLE 
I. 
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TABLE I 

The values of minimum Lorentz factor T, in seconds, E^ ^^^y in MeV used in the 
present work for wind luminosity = 10^^ erg/sec and equipartition parameter 
es = 0.01 from the contour plots given in 0. 



i 


ti^j SGC 




50.12 


0.707 


0.794 


56.23 


0.316 


0.794 


70.79 


1. 


3.162 


79.43 


0.316 


3.162 


89.12 


0.1 


3.162 


1 nn 




7 Q/1Q 

( .y4o 


112.2 


0.316 


19.952 


125.892 


0.316 


125.892 


141.25 


0.251 


125.892 


158.48 


0.63 


79.432 


177.82 


0.251 


79.432 


199.52 


0.501 


50.118 


223.872 


0.125 


50.118 


251.188 


0.1 


31.622 


281.838 


0.398 


19.952 


316.227 


0.079 


19.952 


354.81 


0.63 


7.943 


398.107 


0.199 


7.943 



The mean isotropic equivalent energy released by GRBs is about 10 erg from [|13 
The total wind energy is related to the energy Egrb emitted in neutrinos by the GRB. 
If the energy released in neutrino emission Egrb is 10^^ erg then the wind energy 
is expected to be of the order of lO^^'erg since neutrinos carry about 5% of proton's 
energy and from one proton two muon neutrinos are expected to be produced. 

In Fig.l. of the present work the number of secondary muons from a GRB in km? 
area detector has been shown for different values of the Lorentz factors F and E^ j^^y. 
The energy emitted by the GRB in neutrinos is lO^^erg and the burst is at zenith 
angle 9 = 180°. We have used the value of the threshold energy of IceCube detector 
Eth = lOOGel^. One can calculate the neutrino signals for higher threshold energies 
of the detector. The values of redshift z of the GRB used in the calculation are 
0.05,0.5,1. and 1.5. The Lorentz factor F and the observed photon spectral break 
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energy E^ j^^y together determine the value of the neutrino break energy E^. 

If r is increased the break energy of neutrino spectrum E^ also shifts to higher 
value, if the photon spectral break energy does not increase. The neutrinos produced 
from a GRB of high Lorentz factor and low E!^ j^^y will have their energies higher 
than those produced from a GRB of lower Lorentz factor and higher E^ j^^y. The 
atmospheric background noise decreases for higher threshold energy of the detector. 
So it will be easier for a detector of higher threshold energy to detect these GRB 
neutrinos emitted from GRBs of high Lorentz factors and low photon spectral break 
energies. 

In Fig. 2. of this paper the effect of varying the redshift of the GRB on the production 
of neutrino events in the IccCubc detector has been studied. As stated earlier the 
burst is at zenith angle 9 = 180°, total energy emitted in neutrino production is 
Egrb — lO^^erg, wind luminosity of the GRB is 10^^ erg/sec and threshold of the 
detector is E^h — lOOGeT^. If the distance of the GRB from the observer is large the 
observer receives very few signals. In Fig. 3. of this paper the numbers of neutrino 
induced muons have been calculated for different rcdshifts when the GRB is at zenith 
angle 6 = 92.86°. The rest of the parameters are same as in Fig.2. 

In Fig. 4. of the present paper the dependence of secondary muon events on the energy 
of the GRB Egrb has been shown. If the total energy released by neutrino emission 
Egrb is increased obviously there will be more neutrino generation. The number of 
secondary muons from the GRB increases linearly with Egrb-T^^^ wind luminosity 
of the GRB is assumed to be — lO^^erg/sec and redshift of the GRB is z — 1.5. 

Using the procedure discussed in section 3 we have calculated the atmospheric neu- 
trino generated background muon events per second in IceCube detector for threshold 
energy Eth = lOOGe^. In case of cos^ = —1 we find 0.018 atmospheric neutrino in- 
duced muon events from + per second per km^ area and for cos^ = —0.05 the 
corresponding number is 0.04. 

We denote the duration of time for neutrino emission from a GRB by ti, at the 
source in observer's reference frame. From Fig.2. and Fig. 3. of this paper it is clear 
that a nearby burst of a few seconds duration (t;^) with energy Egrb of the order of 
lO^^erg will be always detectable in IceCube for threshold energy Etu = IQQGeV . 

We have also investigated for a known redshift z and Lorentz factor F of the GRB 
for what range of values of the total energy Egrb the IceCube detector mihEth = 
lOOGeV will be able to detect the GRB at a zenith angle 6. By detectable we mean 
the number of secondary muons produced by the GRB is more than the background 
atmospheric neutrino induced muon events integrated over the observed duration of 
the burst. If the duration of neutrino emission by the GRB is long the background 



10 



noise is more compared to the case when the duration of neutrino emission is short. 
In calculating the observed duration (duration for neutrino emission ) of the burst on 
Earth the correction due to the redshift of the source has been taken into account. 
The observed duration of the GRB from Earth for neutrino emission is x [1 + z) 
where z is the redshift of that GRB. We have calculated number of neutrino events 
from a GRB for cos 6 = —1 and cos^^ = —0.05 and duration of neutrino emission from 
the GRB {ti, ) is assumed to be 1. second. The correction on the matter traversed 
by neutrinos due to 2km depth of the detector below the surface of the Earth has 
been taken into account. In Fig.5. of this paper the regions above the curves are the 
allowed regions of parameters (T and Egrb) for detectable GRBs at redshift ^ = 1.5 
in IceCube detector. One can extend this study for any redshift z, zenith angle of 
observation 9 and duration of neutrino emission {ty) from the GRB. 

Simultaneous observations of photons and neutrinos from a burst in a small time 
window will confirm its signature. The GRBs are distributed uniformly as standard 
candles throughout the Universe |T^. There is a large number of faint GRBs from 



which less than one event will be detected on average in the muon detector. 

From redshift z = 1.5 for F = 100 and Egrb = lO^^erg we can expect 14.4 muon 
events for E^^ = 7.94:3MeV at zenith angle 6 = 180° in IceCube neutrino telescope 
for threshold energy Eth = lOOGel^. Even from a distance of redshift z = 2.1 if total 
energy emitted in GRB neutrinos is lO^^erg we can expect 8 muon events in IceCube 
for threshold energy Efh = lOOGeV at zenith angle 6 = 180° when F = 100 and 

= 7M3MeV. 

So there is possibility of neutrino detection from bursts at high redshifts if they 
emit a large amount of energy in neutrinos Egrb in a short duration of time (t,^) . 



5 Conclusion 

We have investigated on the capability of the next generation kilometer scale IceCube 
detector in detecting high energy neutrinos from individual GRBs. It is always easier 
to detect nearby GRBs. Although the neutrino events steeply fall with increasing 
redshift of the GRB there is possibility of neutrino detection even from the far away 
bursts in IceCube neutrino telescope. If the GRB is of short duration ty, low Lorentz 
factor F, emit large amount of energy Egrb in neutrino emission and if its photon 
spectral break energy E^^j^j^y is high the possibility of detecting neutrinos from that 
GRB increases. Not only understanding the physics of GRBs, testing special relativity 
and equivalence principle but we can also look forward to study physical processes 
occuring in the early Universe with neutrino telescope of km? area. 
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Figure 1: Number of secondary muons in IceCube detector from neutrinos produced 
from a GRB as a function of the Lorentz factor F of the GRB have been plotted in this 
figure. The detector area is Ikrri^ and the GRB is at a zenith angle of ^ = 180°. The 
wind luminosity of the GRB is assumed to be lO^^erg/sec for equipartition parameter 
cb = 0.01. The GRB parameters F, observed variability time ty sec, phton spectral 
break energy E'^^j^^y used here have been given in TABLE I. The total energy released 
in neutrino emissions is assumed to be lO^^erg. 
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Figure 2: The number of neutrino induced muons in km"^ area from a GRB have 
been plotted as a function of redshift z of the GRB. The detector threshold is E^h = 
lOOGel^, the burst is at = 180° and wind luminosity of the GRB is assumed to 
be lO^^erg/sec for equipartition parameter = 0.01. The total energy released in 
neutrinos is assumed to be Egrb = lO^^erg. Here F is the Lorentz factor, is the 
wind variability time and E^ j^^y is the observed photon spectral break energy 0. 
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Figure 3: All the parameters are same as in Fig. 2. only in this case the GRB is 
observed at zenith angle 9 — 92.86° from the South Pole. 
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Figure 4: The dependence of number of neutrino induced muons from a GRB in 
IceCube detector on the total GRB energy Egrb emitted in neutrinos have been 
shown in this figure. The detector area is Ikm'^ and the GRB is at ^ = 180°. The 
redshift of the GRB is assumed to be 2; = 1.5, wind luminosity lO^^erg/sec and 
€3 = 0.01. The threshold energy of the detector is Eth — lOOGeT^. The parameters 
r, ty and E!^ are same as in Fig.2. 
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Figure 5: In this figure for redsliift z = 1.5 and for zenith angles 6 = 180°, 6 = 92.86° 
the minimum total energy emitted in neutrinos by the GRB for which the GRB will 
produce signal above the atmospheric neutrino generated background in km^ area 
detector for minimum Lorentz factor F has been plotted. The threshold energy of 
the detector is assumed to be Efh = lOOGeV. The regions above the curves give the 
allowed regions for positive observations. Here the duration of time (t^) for neutrino 
emission from the GRB is assumed to be 1. second. The wind luminosity is assumed 
to be 10^^ erg/sec and equipartition parameter €b — 0.01 . The values of F, E^^j^^y 
and sec used to do the calculations are given in TABLE I. 
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